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ABSTRACT 


A numerical method for predicting the movement of tropical storms 
has been developed utilizing smoothed 700-mb (Atlantic area) and 500-mb 
(Pacific area) height fields as produced by the Fleet Numerical Weather 
Facility, Monterey, California. Geostrophic steering components are 
computed from the contour analyses at or near forecast time for pre- 
dicting storm positions up to 72 hours. Recent-history vector forecast 
errors are employed as corrections to improve the basic steering fore- 
cast for periods up to 36 hours. Testing of storms from 15 August to 
1 November 1965 in both Atlantic and Pacific regions indicates the 
method is comparable in accuracy to official forecasts as published by 
Fleet Weather Central/Joint Typhoon Warning Center, Guam, and Fleet 
Weather Facility, Jacksonville, Florida, as well as to Atlantic-area 
forecasts made with the NHC-64 statistical approach for periods up to 


48 hours. 
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is Introduction 

Many methods for forecasting tropical storm movement have been 
developed in the past ten years. Presently the most successful are the 
statistical approaches (1).. Some examples are AROWA, Travelers-1966, 
Wang, Fairless, and’ NHC-64 (2)(3). Of increasing accuracy are the numeri- 
cal-models used-at Fleet Numerical Weather Facility (FNWF), Monterey, 
California, and at the Environmental Sciences Services Administration 
(ESSA), Suitland, Maryland. These methods employ parameters such as 
700-mb and 500-mb steering, persistance, climatology, thickness and sea 
level pressure parameters. Most desirable is a dynamic approach to 
‘forecasting tropical storm movement but at present dynamic steering 
models cannot’ fully compete with the statistical and re ess ap=- 
proaches. This is, in part, due to the fact that data networks in 


tropical areas are not ‘dense enough. 


ime term tropical storm *s used throughout the text includes 
depression, storm, and hurricane/typhoon stages. 
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ee Background 

In October 1966 Lt. J. E. Kerr, Naval Postgraduate School (PGS), 
Monterey, California, investigated the feasibility of a geostrophic 
steering model-using FNWF's so-called SR fields (see Section 4) at 
various mandatory~levels’ (4). His method was, in fact, similar to the 
geostrophic steering program used operationally by FNWF (5), but with all 
empirical and climatological factors eliminated. 

Kerr's investigation led to results for 12- and 24-hour forecasts 
which appear competitive with the official forecasts as published by 
Fleet Weather Central, Guam, Mariana Islands (2), and Fleet Weather 
Facility, Jacksonville, Florida (3). 

Specifically, Kerr's work indicates that for 12-hour forecasts the 
Mi o: steering is best for the Atlantic hurricanes and the 500-mb 
SR steering is best for the Pacific typhoons. His 24-hour experiments 
show the 500-mb steering to excel when modifiec by a statistically de- 
rived correction based on direction of movement and position of the 
storm. Kerr's results are presented in tables 1 and 2. All averages 
were computed relative to each storm with no weighting for number of 
forecasts made. However, if the average of the individual forecasts is 
calculated and compared to the official forecast errors, computed simi-~ 
larly, then the accuracy of Kerr's approach is highlighted still further. 
Table 3 displays Kerr's figures alongside the recomputed averages. As 
an example, for 91 le-hour forecasts in the Atlantic, using the SR 1000- 
mb field, Kerr's average vector error becomes 99 nm vice 101 nm while the 
official forecast error for 137 forecasts improves to 86 vice 94 nm. 


Other comparisons are self-evident. 


5. Aims 

In view of Kerr's results and recommendations (4) and the demands 
of operational forecasting the aims of this study may be stated as 
follows: 

a. Search for the optimum geostrophic steering level. Include 
the {00-mb level in such a search in view of the disparity of results 
regarding the best level in the Pacific and Atlantic areas. 

b. Extend the forecast period, using the optimum forecast level, 
to 36, 48 and 72 hours and compare the results to the official forecasts 
as well as to a single outstanding statistical forecast approach (as 
NHC-64 in the Atlantic area). 

c. Find a statistical modification to the results in a and b 
above in order to closely approach a plateau of forecast skill commen- 
surate with data density, accuracy of observations, and methods of nu- 


merically processing input data and computed parameters at FNWF. 


1) 


4. Numerical Procedures 

The SR fields may be viewed as long-wave or macroscale height pat- 
terns which result from a mathematical decomposition, as developed by 
Holl (6). An SR height (2.5) is related to the observed height, Z, as 


follows: 


ye a (1) 


SR —96 


where Lan refers to the disturbance height field which is closely related 
to the short or minor wave components of the height field. For the pur- 
pose of tropical storm forecasting the storm may be considered an SD 
feature, thus the SR field is tantamount to the residual height field 
which results upon removing the storm's circulation. 


The geostrophic steering edquation as applied to the SR field may be 


stated as 


WG = Hp (Mset - ees ) = 


In finite difference form,equation (2) becomes 
—> J44 T-41, A J J \"% 
s- Zenm- P=—(Zee- Z mM 
VG ai) ( SR ZseJA (Zse Z:0)) | as 


An example of using this equation on a linear grid is shown in figure (1). 
Here the grid points, I and J, are equivalent to that in the FNWF opera- 
tional grid but I and J points are not necessarily located at whole- 
numbered FNWF grid points. The storm location to the nearest 1° latitude 
and longitude determines position I, J. 

To avoid division by very small numbers in equatorial latitudes, 
the sine function, as used in the coriolis parameter, f, is modified to 

a[ (sine +L) - Sin 6 | 
+ 45 ¢} 

The graph of this function is shown in figure (2). 
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The steering program was run on the PGS CDC 1604 computer. The 
input data were best track latitude and longitude of the tropical storm 
at 0600 and 1800 GMT (I, J on figure 1) except for Atlantic storm Carol 
for which 0000 and 1200 GMT positions were used. The SR field at 1200 
or 0000 GMT (i.e. six hours after time of storm positions, except for 
Carol) was selected and the four points I, J+l; I+l1, J; I, J-1; I-1l, J 
were located at standard FNWF grid distances from I, J. The meridional 
and zonal geostrophic wind components were then computed by equation (4) 
for each of the four points. The storm (I, J) was steered for one hour 
with the mean of these components. This process was repeated for the 
required number of one-hour time steps, for forecasts at 12-hour intervals, 
up to 72 hours. -Only the single SR analysis selected for the input time 


was utilized. 


ant 


58 Testing and Results 


The procedure outlined in section four was used to forecast the 
movement of the 1965 tropical storms for the Atlantic and Pacific areas 
for intervals up to 72 hours. All storms from mid-August to 1 November 
1965 were processed. However, because of technical difficulties, re- 
sults for all forecast intervals and techniques developed in this study 
are not available for all tropical storms in this period. This period 
includes all storms in the Atlantic area for the year 1965. 

Since Kerr computed forecasts for 12 and 24 hours using SR 500, 
this study first concerned itself with the extension of such forecasts to 
36, 48 and 72 hours (see table 4). The figures are the mean vector 
errors, in nautical miles, with the official mean vector errors listed 
for comparison. Kerr's 12- and 24-hour forecasts are included. The 
figures in paranthesis are the number of forecasts used to compute the 
average. 

The Atlantic tropical storm summary has individual storm official 
errors for 12- and 24-hour forecasts only, while overall averages are 
given for 12-, 2h-, 48- and 72-hour forecasts (3). In the case of 
typhoons, individual and overall average vector errors were published for 
24, 48 and 72 hours (2). 

At this point of the investigation the accuracy of the SR fore- 
casts approach but do not show skill relative to the official forecasts, 
with one exception, the 72-hour Atlantic forecasts. However, the trend 
of accuracy of SR 500 forecast errors relative to official forecasts 
improves with increasing forecast interval, especially for the Atlantic 
storms. 

Since Kerr indicated relative success in the Atlantic using SR 
1000-mb fields for 12-hour forecasts, it was decided to experiment with 
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SR 700, comparing results to 1000-mb in the Atlantic Ror lée-hour fore- 
casts and to SR’ 500 in the Pacific for all time intervals. With 
reference to the 12- and 24-hour forecasts using SR 700 in table 5, it 
became apparent that the 7O00-mb SR fields did not fare as well as the 
500-mb' SR fields:‘in the Pacific, with the ratio of the average le- 

(24-) hour SR 500 to SR JOO 80/100 (172/183) nautical miles. Bess at 
l2hours and Mary at 24 hours provided exceptions to the average result. 
In the Atlantic the SR 700 excelled SR 500 with the ratio of the average 
12-(24-) hour SR 500 to SR 700 108/91 (192/177) with Carol a slight 
exception to the average result for the 2h4-hour forecast. 

Since Kerr regarded the SR 500 to be better than the SR 1000 for 
2h-hour’ forecasts in the Atlantic and in view of the results shown in 
table 5, the 700=mb SR steering was extended to 72 hours in the Atlantic 
only. Results are shown in table 6. 

Upon comparing table.6 to table h, it eee that although SR 700 
is, overall, an improvement over SR 500, the relative accuracy of the 
former to the latter decreases consistently from the le@to the 72-hour 
forecast. In fact by 72 hours the SR 500 is better (489 nm) than SR 700 
(495 nm) but the sample is not homogeneous so the trend may not be valid. 

In attempting to develop an objective method for statistically im- 
proving the SR forecasts, the track derived from: prognostic positions was 
‘compared with the best track. It was apparent that the tracks were 
similar in shape but differed in coincidence of forecast and verifying 
positions (figure 3). This observation suggested application of a cor- 
rection to the basic SR forecast using the recent-history vector forecast 
error: The technique is described symbolically in the following 


equations. 
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aS 
The vector error (E) for a 12-hour forecast as dependent on best 
track (T) and forecast (F) positions may be represented as follows: 
—, es > 
= T 
ay a ay (5) 
where m and z refer to the meridional and zonal components of the error. 


In greater detail, 


— A 
Cong [Te ~T,,) ~ (Fe - Taz demas” (6a) 
and 
ans A 
Ea lo T)- CFo-Taag tan 
Thus, 


| oa =O es oe (7b) 


Zi 


t 
Modified l2-hour forecast storm positions, Fio> starting from ; 
initial conditions at time "O",are represented by equations (8a) and 


(8b). , 


Gratz = Frnagt (To - Fo)miz (8a) 
Fr 29,° Ferg > (Tom Fo) 212 (8b) 


Forecasts for 24, 36, 48 and 72 hours employ the 2h-, 36-, 48- and 72- 
hour recent-history vector error corrections, respectively. See figure 
(4) for a vector diagram of a 12-hour modified forecast. 

Application of such a correction showed considerable improvement in 
the 12-, 24- and 36-hour forecasts. As noted in tables 7, 8 and 9, the 
average errors dramatically decrease from 73 to 46 nm, 172 to 105 nm, and 
248 to 185 nm for 12-, 24-, and 36-hour forecasts in the Pacific, and 
from 91 to 49 nm, 177 to 126 nm, and 262 to 251 mm for 12-, 2h- and 36- 


hour forecasts in the Atlantic. Every storm except Mary at 12-hours 
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showed improvement in the first two forecast intervals. Anna, Betsy, 
Debbie, and Virginia showed deterioration with the recent-history cor- 
rection at 36 hours. At 48 hours the overall average vector error for 
the modified forecast was smaller (table 10) in the Pacific but higher 
in the Atlantic with most storms behaving like the average. At 72 hours 
(table 11) the overall average error with the history-error correction 
did not improve the Atlantic or Pacific forecasts. 

um an effort to improve the longer range forecasts, especially 
where the largest forecast errors occurred during the recurvature period 
the most recent 12-hour errors were computed as in the above equations 
and applied to the 48-hour forecasts in appropriate 12-hour multiples of 
the forecast intervals. In the case of a 48-hour forecast the vector 


error is applied as in equations (9a, 9b) below. 


M? 
Bayes Fryct 4 ( To ~ Fo) maa (9a) 
Fey = Fey, tt CTo- Fo)oan (9b) 


This approach has shown some improvement in reducing the overall vector 
errors in forecast positions, but frequently the tracks connecting prog- 
nostic positions become erratic. Figure 5 shows an example for Carmen. 

Table 10 shows the improvement in the history-error modifications 
by emphasizing the corrections due to the more recent error (F'). 

In the Atlantic the Fo error (345 nm) closely matches the F error 
(346 nm) while in the Pacific an po error of 180 nm is a marked improve- 
ment over F and F with 269 and 230 nm,respectively. 

A logical extension of the above modification technique involves a 
weighting of short-term and long-term vector forecast errors in order to 


incorporate the most recent trend into the forecast along with the longer- 


term corrections. 


ell 


In the case of a 48-hour forecast, application to the meridional 


component may be symbolically stated as follows: 


Fee mvp + 4a (T-F FE nag? Wet(To- 7 + ae(T Fe Fades (10) 


where i and i are weighting factors Rant on si Nocteani interval, 
48 hours, and the intervals 12, 24 and 48 hours, from which the errors 
are computed. As an example, weighting the le-hour forecast error by 

AS. (=W,), the 24-hour forecast error by ys (=W,), and the 48-hour fore- 
cast error by l (=W)), while adjusting the weights for the 48-hour 


~ 


forecast interval by setting W — corer. 


+ I] E J+oeelE, ror |E, wl (11) 


eet 
Limited results using this kiss for a om forecast, Fig » @lve a 


00d 


nye PONS 


mean vector error of 228 nm for 10 forecasts of typhoon Carmen. This 


‘7 
figure does not represent an improvement over F).9 (see table 10). 


ce 


bre Comparison to NHC-64 

In comparing the average vector errors of the éptaiuiant forecast 
method of this investigation to the NHC-64 method, the subject method 
did better than NHC-64 in the 12-, 24- and 36-hour forecasts but worse 
than NHC-64 for 48-hour forecasts. Table 12 shows the results of this 
comparison. Forecasts from hurricane Carol by the subject method and 
NHC-64 were all at 0000 or 1200 GMT. NHC-64 forecasts for other storms 
were also made at 0000 or 1200 GMT while those from the subject method 
were at 0600 or 1800 GMT. For comparison to NHC-64 in the latter case 
the average of the 0600 and 1800 GMT forecasts on either side of the 
NHC-64 forecast was used. Although, as a consequence, the samples are 
not completely homogeneous, the results are similar if Carol alone or 


all storms are considered. 


“optimum forecast method in the Atlantic for this comparison was 
700-mb SR steering with the F' recent-history correction for 12-, 2h- 
and 36-hour forecasts and the unmodified 700-mb SR steering for the 48- 
hour forecasts. 
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kc Concluding Remarks 

The numerical method of tropical storm steering as described here 
gives results which are comparable to or excel official forecasts for 
12-, 24-, 36-, 48- and 72-hour forecasts. Moreover, with reference to 
the most accurate single statistical technique in the Atlantic area, 
NHC-64, the subject method shows excellence for the first 36 hours of 
the forecast interval. It should be nesea that the method just outlined 
is strictly objective employing only analyzed SR fields with the cor- 
rection due to the recent-history error being a substitute for the prog- 
nostic SR field by capturing the most recent trends of storm movement due 
to changes in the SR pattern. Also, use of the history modification 
corrects for use of an improper steering level as there is little doubt 
that a single such level is not appropriate for all storms. 

Two pertinent facts concerning comparison of vector errors should 
be mentioned. One, the forecasts are biased in favor of this investiga- 
tion because the initial input for the numerical computations employed 
the best track position while official forecasts were made from opera- 
tional positions. The average distance of reconnaissance aircraft fixes 
from best track positions was 13 miles for all storms in the 1965 season 
in the Pacific (2) and considerably less in the Atlantic. Two, the 
difference between operational and best track positions may be considered 
of minor significance compared to the fact that the operational "12, 

24, etc." forecasting period is generally longer than the research "12, 
24, etc." forecast period. For example, Fleet Weather Central, Guan, 
files warnings every six hours (00, 06, 12, 18 GMT); the message con- 


tains a warning-time position which really represents a 3- to 12-hour 


ou 


extrapolation of the storm position using the last radar or reconnais- 
sance fix, or most recent surface or upper-air reports. Thus a 24-hour 


forecast is, in effect, a 27- to 36-hour forecast. 


en 


de Recommendations for Further Study 

1. Apply the SR geostrophic steering technique to prognostic 
SR fields. 

2. Check the validity of optimum results from the 1965 tests by 
further tests in the 1966 season. 

3. Evaluate the subject method in the 1967 season under real-time 
operational conditions. Compare results to the present method employed 
by FNWF. 

4, Stratify storm data by region, season, storm intensity, etc., 
and search for regression relations to improve upon the basic geostro- 
phic SR steering (with or without history modification). 

5. Continue tests of history modification by varying the weights 


in equation (11). 
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Hieeril Caties Setficial SR 500 Modified 























Forecasts SR 500 
Anna alT 149 205 
Betsy 130 Lom 146 
Carol noe OW gee it 
Debbie mS5 149 ae): 
Elena 260 353 265 
Average 187 198 iar 
Typhoons 
Lucy 169 216 236 
Mary 107 111 ee) 
Olive 138 129 141 
Rose 5 139 = eS 
Shirley 231 230 182 
Dr ix #30 187 NSS. 
Virginia 289 200 182 
Bess 106 120 105 
Carmen 148 156 12h 
Della we 196 189 
Average 153 168 154 
Table 2 
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Average vector error Average vector error 

















Hurricanes by storm [from (4)] by individual forecasts 
Kerr Official Kerr Official 

12 Hour SR 1000 1yenL || 5] 94[ 5] 99(91)* 86( 137) 

24 Hour SR 500 Mod. 177{ 5] 187{ 5] 166( 84) 171 (2385) 
Typhoons 

12 Hour SR 500 eagles 80( 109). ' 

24 Hour SR 500 Mod. 154{10] 153[10] 154(100) 148(199) 
1 


Numbers in [{ ] represent number of storms. 


“numbers in ( ) represent number of individual forecasts. 


Table 3 
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ie Hour 24 Hour 
Hurricanes SR 700 SR 1500 SR 700 SR 500 
Anna 70( 8) 82( 6) age) 149( 6) 
Betsy 72) 83( 29) eS) 157( 28) 
Carol 127( 31) 183( 29) 161@ 307 
Debbie 54( 9) 89( 10) d12(§ 9) 149( ap 
Elena 1s) 194( 13) 328( 13) 353( 2 
Avon 91( 62) 108( 58) 177( 89) 192( 85) 
Average” 114( 89) 
Typhoons 
Lucy 98( 16) 84( 12) 2 1) 
Mary 55( 8) 28( 6) 106( 7) la (5) 
Olive TTC 8) 61( 6) 156( 8) 1294 5) 
Rose 117( 9) 34( 9) 204( 8) 139{ 8) 
Shirley 1 Sakis) 12a{ 13) 232( 12) ~° 230( 12) 
Tex 108( 15) 97( 14) | 188( 14) 187s) 
Virginia 108( 7) 104( 6) 202( 6) 200(a'B) 
Bess e1( 12) 8&1) 151( 10)  120( 13) 
Carmen 103( 16) on 11) 186( 15) 1566716) 
Della FOO 13} Co er2) 196( 11) 
Average” 100(117) 80(109) 183( 80) SS (eo) 
Average 172(100) 


1 average for only those storms forecasted by both SR 700 and SR 500. 


2 ; ae 
Average for all storms forecasted by method given in column heading. 


Table 5 


Se 


Pann Lcanes 
Anna 
Betsy 
Carol 
Debbie 


Elena 


Average 


he HOU 


TO( 8) 
We 310) 


54( 9) 


170(14) 


91( 62) 


24 Hour 


127( 7) 
141(31) 
183(29) 
ie 1G) 


328(13) 


177( 89) 


Table 6 


33 


36 Hour 


189( 6) 
205(30) 
276(28) 
LT3( 3) 
471(12) 


262( 82) 


48 Hour 


284( 5) 
273(29) 
363(28) 
220( 6) 


595(11) 


346(79) 


i2.heur 


524( 3) 
399( 26) 


261( 5) 
857( 9) 


495(43) 


i] 


Hurricanes SR 700 F F 





12 ie 
Anna 70( 8) 47( 7) 
Betsy 12 45(30) 
Carol 45(29) 
Debbie 54( 9) 28( 8) 
Elena 170( 14) Tals (1.3) 
Avexapen 91( 62) 49(58) 
facaee: 48(87) 








Typhoons SR 500 


Lucy 84( 12) 59( 8) 
Mary 28( 6) Luc 5) 
Olive 61( 6) 
Rose 34( 9) 34( 8) 
Shirley 122 (a4) 46(12) 
Trix 97( 14) 47(12) 
Virginia 104( 6) 4S( 4). 
Bess 84( 14) 49(13) 
Carmen SHA ey 50(16) 
Della 85( 12) 52(10) 
Average™ 73(103) 46 (88) 
Average” 80(109) 


i] 
1 average for only those storms forecasted by both Fo and Flo: 
“ average for all storms forecasted by method given in column heading. 


Table 7 
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Hurricanes 
a 


Fou 
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Poh 





Anna i | 126( 5) 
Betsy 141( 31) Jeli 29) 
Carol 183( 29) 113(26) 
Debbie LI ee 99( 7) 
Elena 328( 13) . 226(11) 
Average 177( 89) 126(78) 
Typhoons SR 500 
Lucy 2NeG 12) 107( 7) 
Mary Pa 6) 95( 3) 
Olive 129( 5) a2 ( 12) 
Rose 139( 8) 58( 6) 
Shirley 230( 12) 163(10) 
Trix 187( 13) 118(10) 
Virginia 200{ 5) 94( 2) 
Bess 120{ > 13) 65.11) 
Carmen 1S 6k re) 89(14) 
Della 196( 11) aa (7) 
Average 172(100) UES 72) 
Table 8 


Hurricanes SR 700 P36 P36 





Anna 189( 6) 269( 3) 
Betsy 205(30) 230(30) 
Carol 276( 28) ee >) 
Debbie Lae) 202( 5) 
Elena 471 (12) —-436( 9) 


Average 262( 84) een 72) 





Typhoons SR 500 


Lucy 130(¢ 5) 
Mary owt 5) 152( 2) 
Olive 
Rose 199( 7) 100( 4) 
Shirley 335 0a) 308( 8) 
Jesio's 275(12) 236( 8) 
Virginia 318( 4) 329( 1) 
Bess 173(12) 137( 9) 
Carmen 2eailoy 145(12) 
Della 287( 9) 142( 5) 
Average” 248 (75) 185(49) 
Average” 180 (54) 





1 werage for only those storms forecasted by both P36 and P36: 


eens for all storms forecasted by method given in column heading. 


Table 9 











Hurricanes SR 700 Fig Fig Fig 
Anna - 284( 5) 432( 1) 5 86ileel ) 
Betsy 273(29) 853(25) 302(25) 
Carol 363( 28) 349(24) 325(24) 
Debbie 220( 6) 205{ 3) 291( 3) 
Elena 595(11) 695( 7) 558. %) 

Average 346(79) 397(60) 345(60) 

Typhoons SR_500 
Lucy 353.0) 267( 6) 

Mary 293( 4) 

Olive 

Rose 254( 6) La .2) Lifaleee ) 
Shirley 425(10) 4i9( 6) 

Trix 348(11) 438( 7) 

Virginia 418(° 3) 

3536 229(11) 211( 7) 183( 7) 
Carmen 307 (14) 264(10) 186(10) 
Della 366( 9) 169(_ 4) 

hee 269( 31) 230(19) 180(19) 

Average 327(71) 293(42) 

Averamee 334(78) 





q I 
1 average for only those storms forecasted by Figs Fig» and Fig: 


i) 
“average for only those storms forecasted by Fig and Figs 


Seat for all storms forecasted by Fy ig: 


Table 10 


Ex 


Hurricanes SR 700 F F 





M fas Te 

Anna 524( 3) 

Betsy 399(26) 551(20) 

Carol 

Debbie 261( 5) 

Elena 8 mn363( 3 
— 517(35) ___ 619(23) 
Average“ 495(43) : 
Typhoons SR 500 

Trix 4o1( 9) none 3) 

Bess 360( 9) 285( 3) 

Carmen 441(12) 528( 6) 
__ Della 550( 6) 4i( 1) 
Average 4h (36) | 471(13) 





? 
amd F_.., 


laverage for~onily~those storms forecasted by both Pre q2 


“average for all storms forecasted by method given in column heading. 


Table 11 


12 Hour 24 Hour 36 Hour 48 Hour 


SR 700 forecasts for 


Hurricane Carol only 39(15) 105 Gi) 175(15) 330(15) 





NHC-64 forecasts for 











Hurricane Carol only TS Ly S) a See) 215@§) 302(25) 
SR 700 forecasts” 40(48) 12% (42) 210(37) 297 (42) 
NHC-64 forecasts? 85(48) 130(42) 227( 37) 267( 42) 





1 Averages for only those forecasts for which data were available for 


both methods. 
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